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A two-particle space correlation function is derived from the single-particle momentum distribution of 
the emission source. A signal of a first order phase transition is obtained from this correlation function if 
density fluctuations are large. 

PACS number: 25.75.GZ, 12.38.Mh, 24.60.-k, 64.60.-i 



The search for phase transitions is a hot topic in 
atomic and nuclear physics. Several different phase 
transitions are observed/or expected to be observed 
like 3 He to 4 He superfluid, Bose-Einstein conden- 
sation, liquid to gas nuclei and hadronic matter to 
Quark Gluon Plasma (QGP). The QGP phase tran- 
sition is predicted by lattice QCD and several QGP 
signals have been suggested M, although no convinc- 
ing evidence Q has yet been presented. The nature 
of the phase transition of QGP depends on the dy- 
namical quarks. A first-order phase transition is ex- 
pected without dynamical quarks. Simulations of lat- 
tice QCD with dynamical quarks have not yet over- 
come the limitations due to finite lattice size. Present 
results indicate a smooth cross-over between phases 
for two light quark flavours (Nf=2) and a first-order 
phase transition for Nj >30. 

Among these QGP signals there are two based on 
the study of particle average transverse momentum 
{{pt)) and on two particle momentum correlations. 
If the dependence of (pt) on the (pseudo) rapid- 
ity density (dn/drj) of produced particles || is stud- 
ied, the phase transition of QGP to hadronic mat- 
ter could be seen as a plateau followed by a second 
rise in {pt)- However this signal is not sensitive to 
the shape of the single-particle momentum distribu- 
tion. When studying the emission source space-time 
extension through Bose-Einstein correlations, it was 
argued that the two particle momentum correlations 
can signal a QGP phase transition j|] through a very 
large apparent source radius in the 'outward' direc- 
tion. 



Recently much attention has been paid to search 
for disoriented chiral condensates (dec) from single- 
particle momentum distributions and two-particle 
momentum correlations Bjorken et al. pointed 
out that large fluctuations in the pion spectra can be 
expected from the dec § which create large variation 
in the source density. 

The fluctuation of the source density near critical 
point was first discussed by Ornstein and ZernikejjJ. 
In this paper we study the two-particle space corre- 
lations of the emission source through density fluctu- 
ations in heavy ion physics. The two-particle space 
correlation function is derived from the single-particle 
momentum distributions. A signal of a first order 
phase transitions sensitive to the shape of the single 
particle momentum distribution is proposed. This 
signal is not only suitable for searching of the QGP 
phase transition but also for the search of the liquid 
to gas phase transition in the excited nucleus. 

Assume a source at rest with volume V in which 
there are N particles distributed in a set of positions 
r*j, i = 1,2,...N. The normalized N particle prob- 
ability distribution is p(fi, r*2, fjv). The particle 
number density at a position r is defined as, 



N 

%=i 



(i) 



Here 5 represents a Dirac delta function. Thus 
p\(f)dr is the number of particles between r - f+df. 
The integration of the source volume will give 



pi(r)df = N. 



(2) 
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The ensemble average of the particle number density by assuming that the two particle space correlation 
pi(f) is assumed to be independent on the position r only depends on the distance between two particles, 
and denoted as p, The statistical average on the left side could be found 

in the following way. The total free energy F is ex- 

p = J dn...[ dr N J2 s(f- fMn,^, ...,? N ) pressed in terms of the free - ener gy densit y 

i=l , 

F = / Fi{r)dr. (11) 



= N J dfi5(f - n)P(fi) 

= ATP(r). (3) For a given temperature, Fi(r) may deviate from its 

equilibrium value F\ due to the density fluctuations. 

The pi(r) is the local particle density which fluctu- Because Fx{r) should be minimum for equilibrium, 

ates around the p and Pl {r) - p reflects the density Let us expand the difference F x (r) - F\ according to 

fluctuation of the source. The volume integral of p Landau's approach^- 
gives the average number of particles (iV) 

/AA m F 1 (^-F- 1 = l -a[ Pl (r)-p] 2 + h(Vp 1 (r)f + ... (12) 

drp(r) = (N). (4) 2 2 



here parameters a and b are independent of the den- 

The two-particle density correlation function can be ., u , , ,-i , , ^ c .. 

* J sity but may depend on the temperature, l^or a finite 

expressed by, g - ze gource ^ e F our i er integral should be replaced by 



p 2 (f,r')= ^ S(f - fi)5(r' - fj). (5) 



the Fourier series. Neglecting higher-order terms in 
eq.(12), 



The ensemble average of P2(r, r') is obtained in a sim- F — F{V = — JZ( a + ^ 2 )|Pi(^)| 2 - (13) 

ilar way as p(r), 

Based on this expression, the probability of having 



p 2 (r, r') = N(N - l)P(r, r'). (6) the fluctuation \pi{k)\ is given by, 



The integral of P2(r,r') is the mean number of two p(\pi(k)\) = A exp(— (a + bk 2 )\pi(k)\ 2 /2VT). (14) 
particle pairs (N(N — 1)) in the volume. 

The density-density correlation is defined as Here A = (f e-( a+bk2 ^ pl ^ 2 / 2VT dpi{k))~ l which im- 



pi(r)pi{r'). From the formulas above it is found plies, 

\pi{k)\ 2 = . (15) 

Pl{r)pi{r') = p5(f- ?) + p 2 (f, ?). (7) a + bk2 

Using this and taking the Fourier inverse in eq.(10), 

One can now study the density-density fluctuations we obtain two par . ticle space correlation function 

from the density-density correlation, 

f 7^ 7r 7^ T ~}\ . JZ ~~)\ 2 /ox P2(r) - p 2 + pd(r) = — — exp[-(-) 1/2 r] (16) 

[pi{r) - p\[pi{r) - p\ = pd(r-r')+p 2 {r,r ')—p . (8) 4nbr b 



The Fourier transformation of p\ (f) — p is ^' ie parameter a must vanish at the phase transition 

point (critical point T = T c ) with zero net-baryon 
Pl (r) -p= (2vr)- 3 J Pl {k) exp(ik ■ r)dk. (9) density, because 

_ d 2 F l (r) 1 dP 

Where pi(k) is the Fourier component of particle a ~ \ Qp 1 ( r y ' T ~ p 1 (r) dp\{r) 

number density fluctuation which is a complex vari- 
able. Because pi(f) - p is real, p\(k) = pi(-k) and Here P is the pressure of the system. For the first- 
because of the normalization pi(f), pi(k)\^ =Q = 0. order phase transition, the two-particle space corre- 
Equations | and | give lation function becomes long-range, 

| Pl (fc)| 2 = Vp + V J e~ fk - r \pW) ~ P 2 W, (10) p- 2 -^- p ^ + p S(f) = ^-, T = T C . (18) 
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In the vicinity of T c , a is non-zero. Thus the two- The average quantum state is as 
particle space correlation function decreases expo- ^ 

nentially. _ _ _ _ _ _ _ _ 

Eq.(16) shows that two-particle space correlation = i^k' ^k* 2 \(h(n) - n)(n(r 2 ) - n)\k^h> \k 2 k' 2 ) 
function contains information of the first order phase _ -> + + ^ 

transition. Now the way of obtaining such informa- = (hk'^k^l 2^ 2^ 2^ a^ag^^ag^ 
tion experimentally will be presented in accordance kl k>1 k ' 2 fc 2 

with the second quantisation method. fit (n)^,-, {f\)fit (f 2 )4>p (r 2 )\hk' ik 2 k' 2 ) . (24) 

Assume that particles with p = hk move freely 

in a static source of volume V at freeze-out. The D2 is not zero although k x = k' 2 and k\ = k 2 . Let 

normalized single-particle wave function is k = k x and k' = k 2 , then the average quantum state 

of eq.(24) will become, 

<W = ^exp(,fc.r). (19) m = ±^^(l± n0 )n i jl i -2H*-*\ (25) 

k k' 

The operators a+ and a % increase and decrease in Here ig thg number of particles in the quantum 
units of the numbers n,; of particles in the various gtate _ The , + , referg tQ the cage of Boge statistics 
quantum states <j> % . The operators become, and ,., to that of Fermi statistics . After the quan . 

turn averaging the D2 must also be averaged in the 
\J/+(f) = ^ (j>j;(r)at, *f/(r) = E^A-(^')^fc (20) statistical sense, i.e. over the equilibrium distribu- 
tion of the particles in the various quantum states. 

which respectively add and remove one particle from Since particles in different qua ntum states are quite 

point r in the system. The operator #+(r)*(r)dV i s independent) the numbe rs n % and ng are averaged 

the operator of the number of particles in the volume independently 

dV. Hence \I/ + \E' can be regarded as an operator n, 



which represents the particle density distribution in (1 ^ n k') n k ~ n fc') n fc - (26) 

s P ace If one now replace the summation in eq.(25) into in- 

~ ~ ,„ tegration, we obtain the following expression for the 

n = m (r)^{r) = ^^ata %l fi % 4> t : / . (21) mean density . density fluctuation; 

k k' 



(n(ri) — n)(h(f 2 ) — n) 

The diagonal terms of the sum (k = k) give the mean 1 f _ -ik-(fi -M c/-> - \ ,? 

. _ . , . ; , = 7 — / n t e %K {Vl r2, 5(ri - r 2 )dk 

density n m the quantum state considered, {2ir) 6 J 

n = <k\}_^ata l :\(f> j :\ \k > (2vr) 6 7 J k ' k 

1 k - -. 1 = nS(r)±\-^ [nte^dkl 2 , (27) 

here f = f\—r 2 and = e ^ e _^/ T is the Bose/Fermi 
Here operator = ata^ and is the number of distribution. 

particles in the quantum state. The operator which After comparison of eq.(27) with eq.(8) the two- 
contributes to the density fluctuation is particle space correlation function is found to be, 

S = EEWf" (2 3 ) pW^-p 2 = Mj^Jn^dk\ 2 - (28) 



n 

k k'^k 



(27 

For a classical (Boltzmann) system, the o^ai" 



The average density-density fluctuation is marked as ata^ so that 
(n(r*i) — n){h(r 2 ) — n). This mean value is calculated k 
in two stages. First of all, the quantum averaging {n{f\) — n){h(r 2 ) — n) 

over all the quantum states of 4>% and 4>g then the 1 . f _ i%.?,?\ 2 , 0Q ^ 

ensember (statistics) averaging should be carried out. (27r) 6 J k 
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and the two particle space correlation function is 

1 



P2(ri,r 2 ) - p + pS(n - f 2 ) 



(2tt)6 1 J n k elk ' r ^\ 2 - 
(30) 

Here = exp((// — e)/T) is the Boltzmann distribu- 
tion. 

If the spin S of the particles is taken into account, 
the two particle space correlation function will be- 
come 



P2{n,r 2 ) - p 
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(2 %ll J V lk ' T dk\ (31) 

Here a is the spin projection and g = 2S + 1 is the 
degenerace of the particle. ri? is the mean number 
of particles at a state of given spin projection a. 

We must of course be aware of that the difference 
in momentum between two Bosons tend to be small 
and more pairs are observed than for Fermions. 

Equation (28) also shows that the difference be- 
tween Bosons and Fermions is similar in configura- 
tion space, the presence of one Boson at some point 
increases the probability that another Boson is close 
to that point. 

Comparing the experimental distribution of two- 
particle space correlation calculated from eq. 28 or 



30| with that from eq. we could probe whether 
the first order phase transition has taken place. If no 
clasical density fluctuations and no quantum correla- 
tions, the particles will be emitted independently in 
space. Otherwise the correlation function becomes 
long range if particles are emitted at the critical 
point. 

Figure 1 shows the dependence of the two particle 
space correlation on the distance between two parti- 
cles (r) with and without a first order phase transi- 
tion. Curve 1 represents the two-particle space cor- 
relation including the phase transistion (a = 0) and 
curve 2 without the phase transistion (a/6 ~ 4). The 
temperature T=160 MeV is chosen. The calculations 
show that if parameter a is of the same magnitude 
as 6, the difference of the two-particle emission pat- 
tern from that with the phase transition is clearly 
visible. Two-particle emission with small separation 
in distance dominates if there is no phase transition. 
The results are also calculated with a pure Boltz- 
mann distribution with T =160 MeV and Bose dis- 
tributions with two different temperatures. These 
results show no signal of a phase transition as ex- 
pected. A possible single-particle momentum distri- 
bution (strong enhancement at low momentum and 



an enhancement at high momentum) which results 
in a long range correlation of two particles in space 
is also indicated in the figure. Experimentally, one 
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Figure 1: The dependence of two-particle space 
correlation function on the separation distance 
between two particles. 

should search for the phase transition on an event- 
by-event basis, looking for significant fluctuations in 
a single-particle momentum distribution. Assuming 
that there is 30% fluctuation of the magnitude in 
the fitted single-particle momentum distribution at 
momentum less than 30 MeV/c due to detector res- 
olution but keeping the function of the momentum 
distribution the same as that indicated by the solid 
triangle in the figure, the two-particle space correla- 
tion function will vary and the varied region is also 
indicated in the figure. 

Figure 2 shows the single-particle momentum dis- 
tributions which are used to calculate the two- 
particle space correlation. It demonstrates that what 
kind of shape of single-particle momentum distribu- 
tion possibly indicate a phase transition. 



i o 




Figure 2: Comparison of two kind of single- 
particle momentum distributions used in the cal- 
culation of two-particle space correlation. 

In eq.(26) the assumption of independent particles 
in different quantum states simplified the calculation 
of the two-particle space correlation. In fact this 
approximation can be removed and the two-particle 
space correlation can be deduced from single-particle 
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and two-particle momentum distributions. 
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